Translocation time of periodically forced polymer chains 
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We show the presence of both a minimum and clear oscillations in the frequency dependence of 
the translocation time of a polymer described as a unidimensional Rouse chain driven by a spatially 
localized oscillating linear potential. The observed oscillations of the mean translocation time arise 
from the synchronization between the very mean translocation time and the period of the external 
force. We have checked the robustness of the frequency value for the minimum translocation time 
by changing the damping parameter, finding a very simple relationship between this frequency and 
the correspondent translocation time. The translocation time as a function of the polymer length 
has been also evaluated, finding a precise L 2 scaling. Furthermore, the role played by the thermal 
fluctuations described as a Gaussian uncorrelated noise has been also investigated, and the analogies 
with the resonant activation phenomenon are commented. 

PACS numbers: 05.40.-a, 87.15.A-, 87.10.-e, 36.20.-r 



I. INTRODUCTION 



The transport features of molecules and polymers at- 
tract nowadays more and more interest. In particular, 
the translocation features of DNA polymers through pore 
cells have been studied with great attention [l| and dif- 
ferent models have been introduced to describe such pro- 
cess. In between them we can mention for instance the 
single barrier potential and the flashing ratchet mod- 
els Q. Despite all these efforts, the basic understanding 
of the translocation process of polymers is still missing. 
In many cases, the translocation phenomena involve also 
molecular motors 0, Q whose complex action has been 
only recently addressed with high attention. In addition, 
nanotcchnological machines can learn from the biological 
processes and thus to work in a simpler and more effi- 
cient way [|| . The goal of this article is to study a cyclic 
time-dependent model that can depict a simple nanotech- 
nological machine which is able to drive with a sinusoidal 
force a polymer chain in one direction. This simple ob- 
ject constitutes, in our purpose, a basic model that will 
be further developed into a more realistic machine able 
to emulate the basic behavior of a biological motor. We 
find a non trivial behavior of the polymer translocation 
time as a function of the driving frequency v with the 
occurrence of a series of oscillations. These oscillations 
arc not related to the vibrational internal modes of the 
chain but to the driving frequency. The first minimum 
occurs at mean times T m < T m , where T m = l/v m is the 
period of the driving at this minimum. In this sense, the 



translocation time dependence presents some similarities 
with the resonant activation phenomena Q- The poly- 
mer scaling properties have been checked and verified as 
a function of both the number of polymer beads, and 
the damping. The article is structured in the following 
way: next section presents the model, then we report the 
results in Sec. III. In Subsec. IIIc the dependence of the 
stall force of the system is also studied, finding a strong 
nonmonotonic behavior with the frequency even in the 
presence of relatively high damping. We finish with a 
conclusions section. 



II. THE MODEL 
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FIG. 1: (Color online) Scheme of the oscillating pushing force 
acting on a polymer chain formed by N monomers. The force 
drives only a small fraction of the whole chain. 
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The polymer is modeled as a unidimensional chain of N 
dimensionlcss monomers connected by harmonic springs 



2 



|, whose total potential energy is 



k N ~ 1 

Vhar = T X] 0^+1 



do) 



(1) 



Here fc is the chain elastic constant, Xi the position of the 
i-th monomer and do the equilibrium distance between 
adjacent monomers. 

The translocation is helped by the presence of a linear 
potential (see Fig. [1} oscillating with frequency v which 
acts in a space region of length Lm- In particular: 



V sp (x,t) 



-F x [1 - cos(27ri4 + 4>)] , x e [-L M , 0] 
0, otherwise 



(2) 

The dynamics of the i th monomer of the chain is then 
described by the following underdamped Langevin equa- 
tion: 



dV h ar dV s , 



dxi dxi 



The polymer is subject to a viscosity parameter 7, being 
neglected other more complex hydrodynamical effects. 
Thermal fluctuations are described as Gaussian uncorre- 
cted noise & (t) which satisfies the usual statistical prop- 
erties (&(*)) = and (€i(t)£j(t + T)) = 2jDS ij S(t) with 
= 1...N). Fq and v represents respectively the am- 
plitude and the frequency of the forcing field, and (f> its 
initial phase. 

The model here described represents a simple example 
of a driven polymer translocation, which may reveal as 
proper mechanism in nanotechnological devices and as a 
reference example in the area of the translocation driving 
machines. With this we mean a system which is able to 
push a polymer chain from one point to another, in anal- 
ogy with a molecular motor which push a long molecule 
in the biological translocation process. 



III. RESULTS 

In this manuscript we will focus on studying the depen- 
dence of the mean first passage time, as defined below, 
and of the stall force of the polymer with the frequency 
of the external driving. Simulations were done at differ- 
ent values of damping and noise (usually 7 = 1.0 and 
D = 0.01) for F Q = 0.1, k = 1 and L M = 5.5. We have 
studied chains with iV=12, 18,24,36 and 120 monomers. 
The computed mean values are the result of averaging 
over 10 4 realizations. For each case we have integrated 
the Langevin equation of motion of the system by using 
a stochastic Runge-Kutta algorithm [l5| with dt = 0.01, 
which is small enough for all the calculations here per- 
formed. 

The initial spatial configuration of the polymer is put 
with all the beads lying at the rest distance do = 1 with 
respect to each other, and the last sited at coordinate 



xjv(0) = 0, the final position of the external force region. 
Every simulation stops when the coordinate of the center 
of mass Xcm = V^Sf 2 -* °f the chain reaches the po- 
sition x cm = and the mean first passage time (MFPT) 
is then computed as 



1 f 27r 

MFPT = T =—J (r(0))d0, 



(3) 



where r(</>) represent the FPT of a single realization, and 
</> is randomly chosen. 
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FIG. 2: (Color online) Log plot of the mean first passage time 
t as a function of frequency of the oscillating force for the 
damping parameter 7 = 1. The minimum and the oscillating 
behavior is showed in the inset. The thermal noise intensity 
is D = 0.01. The first minimum, which is also the global one, 
satisfies the very robust condition: = T m = ctr m , were a 
is a constant. 

The adoption of the center of mass as reference point 
in the movement has the role of avoiding the boundary 
effects given by the spatial extension of the motor. In 
fact at the end of the translocation process of the full 
chain, the total force acting on the polymer is smaller 
because of the decreasing number of monomers inside the 
motor. The choice of the center of mass as a reference 
will preserve the calculations from this unwanted bound- 
ary effect, being the average number of monomers inside 
the motor always the highest possible. Fig. [5] shows the 
MFPT t as a function of the frequency of the oscillating 
force for damping parameter 7 = 1, which corresponds 
to a quite damped condition for the system. A clear min- 
imum is visible at a resonant frequency v m ~ 6 • 10~ 3 . 
Moreover, a very strong oscillating behavior appears in 
the region v 6 [4 • 10 -3 — 4 ■ 10~ 2 ], before the high fre- 
quency saturation value is reached. The thermal noise 
intensity is given by D = 0.01. Temperature acts here in 
destructive way, hiding the oscillating information con- 
tained into the MFPT output, as we will study later on 
(see in Fig. [8]). 

Fig. El shows the standard deviation a (left panel), 
the mean velocity (center), and the mean elongation 
dmed of the monomers with respect the rest distance 
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FIG. 3: (Color online) Standard deviation a (left), mean ve- 
locity (center), and mean difference between the rest distance 
do and the distance of the monomers during the dynamics 
(right) related to Fig[2] The negative values of d me d means 
that the springs are compressed. The oscillation appears in 
both the mean velocity and the elongation difference. 
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do during the dynamics (right). d mec t is calculated as 
the mean along the chain of the distance between sub- 
sequent monomers, then further averaged in the number 
experiments. The negative values of d me d means that 
the springs are in average compressed. The oscillation 
appears both in the mean velocity and in the elongation. 
The minimum passage time is obtained for the maximum 
velocity which also corresponds to the minimum elonga- 
tion of the springs, i.e. the maximum average compres- 
sion. 

The computed value of r has been averaged over the 
initial phase <f>. In the low frequency limit we can com- 
pute the MFPT value as 



TM 



1 

2^ 



T s {(j>)dcj}. 



Here t s is the value of the MFPT in the case of fixed po- 
tential, being the low frequency limit an asymptotic limit 
where the time scale of the dynamics is much faster than 
the period T = 1/v of the external force. In that case, 
the polymer " sees" a static potential in all his transloca- 
tion time. Because the sliding is very slow when the force 
is close to zero (4> = 0, free diffusion case), the MFPT 
tends to grow up and this explains its high values. 

In opposite way, the high frequency limit corresponds 
to a dynamical translocation scale much slower than the 
period T. In that case the polymer sees a constant mean 
force Fq, being the cosine contribution averaged to zero 
during the slow dynamics. In this case the MFPT can 
be easily evaluated analytically: 



i(*) 



n 7 F 
Nj 



1 



(5) 



FIG. 4: (Color online) Mean first passage time r as a function 
of frequency of the oscillating force for the damping parameter 
7 = 1 with various values of the initial phase 4>, namely <j> = 
0.01, 7r/4, 7r/2, 37r/4, 57r/4, Sn/2, 7n/4. A minimum is always 
present. Moreover a very strong oscillating behavior appears 
in the region before the high frequency saturation value. The 
thermal noise intensity is D — 0.01. 



scendent, and the number n 7 cannot be guessed in ad- 
vance. In the overdamped case (7 — > 00) with also the 
related time rescaling (r' = r/7), the equation becomes: 



L C m( T ) 



n 7 F , 

T 

N 



(7) 



(4) A rough estimation of r' can be obtained using the rest 



value for the number of monomers into the field, n 7 = 5.5 
and L cm as the distance between the initial position of 
the polymer center of mass (x cm (0) = —L crn ) and the 
position x = of the boundary. Remembering the initial 
condition, L cm = (N — l)/2 = 5.5, we have, finally: 



, N(N-l) 

t' = — i '- ~ 120, 

2n 7J Fb 



(8) 



which is quite close to the high frequency value of Fig. [51 
where the overdamped limit is plotted together with 
other curves. 

The intermediate region is characterized by a translo- 
cation time of the same order of magnitude than the pe- 
riod of the external force, r ~ T. In this special case 
the polymer can, in average, reach the boundary dur- 
ing the first period of the driving. Thus the interplay 
between the sliding down time and the external force 
period reaches its maximum synchronization at this frc- 



where n 7 is the mean number of monomers at a given 7 ^ nc y v ™- Th e rerond mitumum is reached when the 
inside the force region at the high frequency limit. Inte- 
grating Eq. [5] in time we get: 



Nj 



1 

1 



yr _ -g 



(6) 



which is the distance covered by the polymer center of 
mass in the time r. Unfortunately the equation is tran- 



polymcr escapes during the 2 nd period of the driving and 
so on. This explains why the effect is a classical and not 
noise-induced: the synchronization with the force driving 
occurs when the boundary is reached at the proper time, 
and this can also occur without any fluctuations, though 
the latter are unavoidable in the dynamics. 

For a deeper understanding of the basic dynamics, a 
set of simulations have been done by fixing the initial 
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FIG. 5: (Color online) Resonant MFPT (r m ) vs resonant pe- 
riod (T m ) for all the case investigated. The data concerns 
the parameters 7 = 0.01,0.1,1,10,15,20 and also the over- 
damped case with ./V = 12. Furthermore, for 7 = 1, the 
different polymer length N = 12, 18, 24, 36. All the points lie 
on the linear relation given by Eq. [9] 



phase, reported in Fig. @] We can see that the MFPT 
can show as well a nonmonotonic behavior, depending on 
the <f> value. According with the previous comments, for 
low frequency values, we observe very high r values for 
4> = 0.01 (very small force during all the dynamics), while 
for high frequency values again the static mean force Fq 
contributes to the dynamic, giving the same r value for 
all the different initial phases, which is also the value of 
the randomly chosen initial conditions. The intermediate 
region shows first minima or first maxima at different 
positions, depending on increasing or decreasing initial 
forces respectively, followed again by oscillations. It is 
worth to note that the couple of opposite initial phases 
{4> = 7r/4 and <fi = 57r/4, or (j> = 37r/4 and (j> — 7ir/4) show 
the same r value in the low frequency limit, because the 
force is the same nevertheless the direction (upward or 
downward) of the cosine function is, and for the couple 
(j) = 7r/2 and <j> = 3ir/2 this limit is the same of the high 
frequency one, because the limit force is Fo in both cases. 

As we mentioned previously, r m , the first minimum 
of the MFPT versus frequency curve, which is also the 
global minimum, satisfies a very robust relation with the 
period of the driving for this minimum, T m , given by 



(9) 



were a = 1.36 ± 0.02 (see Fig. [5]). This means that the 
minimum translocation time is obtained for a driving pe- 
riod which approximately corresponds to one and one 
third the time the polymer takes to cross the boundary. 
The above condition has been checked for all the values 
of 7 and all the polymer length we have studied. 

As a further confirmation of Eq. |H1 we show in Fig. [5] 
the oscillating MFPT region for different values of damp- 
ing. In order to compare the curves, the MFPT and the 
frequency have been scaled by damping. The presence of 
a minimum followed by a series of oscillations is observed 



for all the values of the damping. In fact, by diminishing 
the damping we obtain higher r/7 ratios but the position 
of the minima and the amplitude of the oscillations do 
not change. It shows that the robustness of the oscillat- 
ing pattern goes beyond Eq. [9] 

Further, the presence of the fluctuations with different 
damping at the same amplitudes (in the scaled values), 
demonstrate once more that the oscillating behavior has 
nothing to do with inertial features and inner vibrational 
chain modes, as one can infer at a first look of the affect. 




v * y 

FIG. 6: (Color online) The minimum and the oscillating be- 
havior is present independently on the damping parameter 7. 
The figure show the calculations for different 7 values: 0.5, 1, 
10, 20, plotted following the scaling indicated in the legend of 
the axis. 7 = 10 approaches very well the overdamped limit 
(also shown). In the inset we show the curves without scaling. 



A. Length scaling 

Based on the hypothesis that in a polymer the friction 
is proportional to L M , for the driving case it has been pro- 
posed 0] that the translocation time of a polymer follows 
a scaling law with respect to its length L as r cx L 1+A \ 
where \x is a certain constant. This doubtful dependence 
has been recently discussed in Ref. [uj LIU ■ ^ n our unidi- 
mensional case the scaling law we recover appears to be 
in good agreement with the Rouse prediction d, [I2J 

txL 2 . (10) 

Fig. [7] shows that the scaling properties hold here in a 
precise way. There we show our numeric results for four 
different polymer size, N = 12, 18, 24, 36, 120. It is ob- 
served that the translocation times follow perfectly the 
scaling law r oc (N — l) 2 where N — 1 is proportional 
to the total length: L = do(N — 1) with do the rest 
distance between two subsequent monomers. This rela- 
tion confirms again the robustness of Eq. [9] The scaling 
law given in Eq. [TU] has been also proved by other works 
where a constant force was used to drive the polymer 
translocation in a Id case fiol [l3|. 
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B. Thermal noise contribution 

The oscillating behavior depends on the thermal noise 
in a destructive way. As we can see in Fig. [8j the role of 
the thermal noise then appears to be to destroy the oscil- 
lations. These are present in the classical dynamics (no 
noise), but the number of visible maxima decreases by 
increasing the noise intensity, and the value of its max- 
ima decreases monotonically. Only the first maximum 
remains barely visible in the plot at the highest temper- 
ature values there shown. It has to be observed that the 
behavior shown in Fig. [8] for high values of D, is qualita- 
tively the same of the typical resonant activation (RA) 
phenomenon Q, where a particle overcomes a potential 
barrier. RA effect has been observed also in polymer dy- 
namics |14| with a characteristic curve very similar to the 



behavior shown here. In this sense an experimental curve 
showing similar behavior as a function of the frequency 
to the one here described, can be related to the presence 
of a fluctuating barrier or, as in the case here presented, 
to an oscillating linear driving in an high temperature 
environment. 



C. Stall Force 

A last result of the model here studied, useful for fu- 
ture analysis and comparisons with other models or ex- 
perimental outcomes, is the evaluation of the stall force, 
i.e. the force necessary to stop the polymer transloca- 
tion. In order to do that, a set of simulation have been 
performed by applying a pulling force F p (See Fig. [5]) on 
the left side of the chain. The initial condition has been 
fixed now with the polymer center of mass in the center 
of the external potential, and the velocity of the center 
of mass is measured waiting for the exit on the left or on 
the right of the potential region. A zero mean velocity 




FIG. 8: (Color online) Range of the oscillating effect for dif- 
ferent thermal noise intensity, namely: D — 0,0.1,0.5,0.1. 
The oscillating behavior tends to disappear by increasing the 
temperature. The maxima decrease and the minima increase. 
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FIG. 9: (Color online) A force pull is applied at the first 
monomer to measure the stall force of the potential. 



allows to estimate the stall force of the system. 

As it is visible in the inset of Fig. [5J the mean veloc- 
ity of the polymer shows a strong oscillating behavior, 
similar for all the 7 values here used. 

Surprisingly, the oscillations are not always present in 
the velocity used for the stall force evaluation. In fact for 
7 = 0.1 they are present, while for 7=1 they are not. 
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FIG. 7: (Color online) MFPT as a function of v for different 
values of the polymer length L, namely N = 12, 18, 24, 36, 
and 120. When axis are properly scaled all the original curves 
(shown in the inset) superimpose upon each other. 
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FIG. 10: (Color online) Stall force as a function of the fre- 
quency of the oscillating driving for 7 = 0.1. The other pa- 
rameters are the same of Fig. [21 The inset shows the mean 
velocity versus frequency curves for two values of the pull 
force, F„ = -0.446 and -0.452. 
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The inset of the two figures [TOl and [TT1 show the two cases. 
Similarly the two related stall forces present different be- 
havior. With two minima the one with 7 = 0.1 (Fig. ITU]) , 
with only one the other with 7=1 (Fig. ITT]) . In both 




FIG. 11: (Color online) Stall force as a function of the fre- 
quency of the oscillating driving for 7=1. The two minima 
displayed in Fig. [10] have disappeared, and now only a mini- 
mum appears. The other parameters are the same of Fig. [2] 
The inset shows the mean velocity versus frequency curves 
for three values of the pull force, F p — —0.446, —0.450, and 
-0.454. 

cases a strong mean frequency dependence of the stall 
force is depicted, though the intensity of force variation 
appears relatively not very high. In fact the variation in 
the scale is somehow small, and an experimental verifica- 
tion could be not immediately simple to perform. In this 
sense the more realistic model depicted in [i~o| is much 
more checkable. 



mer driven by a uniform force oscillating sinusoidally 
in time and acting on a limited spatial region. The 
translocation time, as well as the stall force of the sys- 
tem show a strong dependence on the frequency of the 
driving. This dependence is visually similar, in many 
respects, and with many distinguo, to the resonant acti- 
vation phenomenon. The driving frequency for the min- 
imum translocation time follows with strict proportion- 
ality the computed value of this very time. A simple 
scaling law for the polymer length, r oc L 2 has been 
found. The role of the noise intensity has been also in- 
vestigated, which gives, in this case, destructive effects 
on the registered MFPT oscillating behavior. Finally, we 
have studied the stall force of the motor as a function 
of the driving frequency, finding a clear non monotonic 
behavior. 



We present here results for a one-dimensional model. 
Further investigations will consider dimensionality effects 
on the quantities here calculated. Dimensionality could 
be irrelevant for polymers constrained to move in con- 
fined channels. Moreover, in many experimental situa- 
tion the polymer is stretched, thus removing the dimen- 
sionality dependence of the measured quantities. 



The model can have artificial application in nanotcch- 
nological devices driven by oscillating fields, and repre- 
sents a first step to further analysis and applications to 
realistic biological translocation processes. 



IV. CONCLUSIONS 

We have studied the main characteristics of a simple 
toy model for studying the translocation time of a poly- 
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